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We propose and analyze theoretically a cavity optomechanical analog of a heat pump that uses a 
polariton fluid to cool mechanical modes coupled to a single pre-cooled phonon mode via external 
modulation of the substrate of the mechanical resonator. This approach permits to cool phonon 
modes of arbitrary frequencies not limited by the cavity-optical field detuning deep into the quantum 
regime from room temperature. 


Quantum optomechanics has witnessed rapid progress 
in recent years, resulting in the successful cooling of me¬ 
chanical modes to their ground state of motion with a 
combination of cryogenic techniques and o^ical sideband 
cooling realized in a number of systems [ll-lSj . This opens 
up a broad spectrum of applications ranging from force 
and field sensing to tests of the foundations of physics, in¬ 
cluding the exploration of the elusive boundary between 
the quantum and the classical world, studies in quantum 
thermodynamics, and more. 


The recent demonstration of parametric coupling of 
mechanical modes deep in the quantum regime offers the 
opportunity to explore a number of aspects of multimode 
phononics, such as the generation of nonclassical states of 
mechanical motion and the development of phonon inter¬ 
ferometry 0,0. While optomechanical sideband cooling 
00 is typically applied to a single mechanical mode, in 
such multimode applications, and more generally in the 
emerging area of nonlinear phononics [s - [l3 1 there is much 
interest in simultaneously cooling two or more modes of 
relatively arbitrary frequencies |l4l |. 


We show that this can be achieved in a cooling cycle 
that relies on the properties of the normal modes (polari- 
tons) of optomechanically coupled optical and mechani¬ 
cal fields [l5l [13 ■ Depending on the frequency detuning 
between the optical and mechanical fields these are re- 
doniinantly either photonic or phononic. In the photonic 
limit they are coupled to a thermal reservoir that is es¬ 
sentially at zero temperature, and in the mostly phononic 
state to a reservoir at the temperature of the mechani¬ 
cal substrate. In that limit they can also be easily me¬ 
chanically coupled to any other mode of oscillation of 
the mechanics by forcing an external modulation of the 
substrate of the mechanical resonator, as was experimen¬ 
tally realized in Ref. 0 |. The advantage of this approach 
is that specific requirements for the optomechanical cou¬ 
pling and frequencies apply only to the two modes (op¬ 
tical and mechanical) comprising the polaritons, with no 
restrictions on the additional mechanical modes to be 
cooled. 


The proposed cooling cycle uses a precooled polariton 
mode as a “polariton fluid” whose nature is first changed 
from photon-like to phonon-like by controlling the pump- 


cavity detuning. When in the phonon-like state it is 
parametrically coupled to the phonon mode to be cooled 
for a duration such that an approximate coherent state 
transfer is achieved between them. The detuning is then 
adiabatically returned to a value for which the polariton 
is photon-like. Its thermalization at the temperature of 
the photon bath (with « 0 and thermal occupation 
ria ~ 0 at optical frequencies 0 ) irreversibly dumps the 
excitations that it carried away from the mechanics to 
the environment, thereby completing the extraction of 
energy from that mode. 

This system is reminiscent of a non-equilibrium heat 
pump |l8| . with the important difference that instead of 
changing the temperature of the cooling fluid, we mod- 
ify its environment via the polariton dispersion relation 
[19| . The role of the expansion phase of the refrigerant is 
achieved by changing the polariton fluid from photon-like 
to phonon-like, and the compression-like phase by the re¬ 
verse process . The heat exchange between the mode to 
be cooled and the fluid is achieved by phonon population 
transfer, and heat disposal is achieved by cavity dissipa¬ 
tion of the polariton fluid in its photon-like form. 

Model. We consider a generic arrangement where a 
single mode of an electromagnetic cavity is coupled op¬ 
tomechanically to a single vibration mode of a mechani¬ 
cal resonator, which in turn can be coupled to a second 
mechanical oscillation mode via external actuation, with 
total Hamiltonian 


H = Horn + + Hdi, (1) 

where i?om accounts for the radiation pressure coupling 
between the cavity mode and one of the mechanical 
modes, Hst describes the interaction between the me¬ 
chanical modes, and iJdis accounts for the intracavity 
field and mechanical oscillators damping of rates k and 
7, with K 3> 7 typically. 

The optomechanical component of the system consists 
of a cavity mode of frequency Wq coupled to a mechanical 
mode of frequency LOb- It is driven by an optical held of 
frequency ujp and amplitude ain. We assume that the 
intracavity held is strong enough that it can be described 
as the sum of a large mean held a and small quantum 
huctuations. In a frame rotating at Wp the Hamiltonian 
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-ffom can then be linearized as 

iJom = —hA{t)a'^d + huJbb^b + hg{b + b^){d + d^) ( 2 ) 

where the bosonic annihilation operators d and b ac¬ 
count for the fluctuations of the optical and mechani¬ 
cal mode around their mean amplitudes a and (3. The 
radiation pressure interaction is quantified by the con¬ 
stant g = ago, with g^ the single-photon optomechan¬ 
ical coupling strength, and the laser-cavity detuning 
A(f) = (ujp — uJa)(t) — 2/3( 7 o includes the mean radiation- 
pressure-induced change in resonator length. In steady 
state and for small damping, we have a w ain/A and 
/3 ~ —g^a^juJb- We take a to be real, which can be 
achieved for an appropriate phase of the pump field. 

The coupling between the phonon mode ‘b’ at fre¬ 
quency (jJb and a second vibration mode ‘c’ at frequency 
uic can be realized by actuating the substrate at a fre¬ 
quency iOs close to their frequency difference, Wg ~ 
LOc — uJb 0- The resulting coupling is described in a 
rotating frame at the modulation frequency ojg and in 
the rotating wave approximation by a "state transfer" 
Hamiltonian 

Hst = hSc^c -I- -I- c^b) (3) 

where c is the annihilation operator for mode ‘c’, S = 
ujc — ujs — uJbi and fio is the parametric coupling strength 
of the two modes proportional to the amplitude of oscil¬ 
lations of the substrate. 

In the absence of mechanical coupling, Hg = 0, the sys¬ 
tem is driven by the familiar linear optomechanical inter¬ 
action only. We focus on the red-detuned regime A < 0, 
which in general leads to stable dynamics for small de¬ 
cay rates, and perform a Bogoliubov transformation to 
diagonalize idom in terms of polaritons described by the 
bosonic annihilation operators A and B [I1[I1|. This 
gives 

Hom = buvA(A)AU + hu;B(A)B^B, (4) 

with frequencies 

P ^_-I 1/2 

^ A ^ A2 -h ± \/(A2 - a;^)2 - Wg^AuJb 
uja,b{A) = - - - 

A plot of uja,b(A) can be found for instance in Ref. |15| . 

The photon-phonon polaritons are coherent superpo¬ 
sitions of the cavity field ‘a’ and the mechanical mode 
‘b’ [^. For A <C —u!b, the low-energy polariton branch 
‘B’, characterized by the bosonic annihilation operator 
B and the frequency ojb(A), describes phonon-like exci¬ 
tations, with ojb approaching Wb- In contrast, on the 
other side of the avoided crossing, —ujb A < 0, 
and in the weak coupling regime g/ujb 1 , the oper¬ 
ator B annihilates photon-like excitations of frequency 
wb —A. The opposite holds for the polariton branch 


‘A’, which is photon-like for frequencies far red-detuned 
from A = —Wb, and phonon-like near cavity resonance. 

Cooling cycle. We now explain how the polariton 
branch ‘A’ can be exploited as a quantum heat pump 
refrigerant to cool an arbitrary mechanical mode ‘c’ de¬ 
tuned in frequency from the bare mechanical mode ‘b’. 

We assume that the photon and phonon modes ‘a’ 
and ‘b’ are initially thermalized at temperatures ss 0 
and Tb [2l|. The cooling cycle comprises four steps: (1) 
adiabatic change in the frequency loa of the polariton 
fluid, with an effect similar to the expansion step in con¬ 
ventional heat pumps, and loosely called adiabatic ‘ex¬ 
pansion’ in the following for that reason; (2) “heat ex¬ 
change” between the flnid and the mechanical mode to 
be cooled; (3) adiabatic change in uja to achieve an ana¬ 
log of the compression stage (adiabatic ‘compression’); 
and (4) thermalization of the fluid to Ta ~ 0. 

The ‘expansion’ is realized by adiabatically changing 
the detuning A{t) from a large negative value A^ —Wf,, 
where the polariton ‘A’ is photon-like (A r::; a), to a small 
negative value —uib <C Aj < 0 close to 0, where it is 
phonon-like {A Ki b). The mechanical coupling is absent 
at this stage, Hg = 0. This transformation conserves the 
initial thermal populations of the two polariton modes: 
the polariton ‘A’ remains unpopulated since « 0 while 
the polariton ‘B’ stays at the temperature of mode ‘b’. 
Since the polaritons have then exchanged their nature, 
this transformation effectively swaps the thermal pop¬ 
ulation between photon and phonon fluctuations. The 
duration n of this step should be slow enough to guar¬ 
antee the adiabaticity of the evolution, ri ^ 1 /( 25 ), but 
fast compared to the photon and phonon damping times, 
Tl < 1 /k, 1/7. 

Once the detuning has reached the value A/ the cou¬ 
pling Hg(f) between the mechanical modes ‘b’ and ‘c’ is 
switched on for a duration T 2 . (For simplicity we consider 
a square pulse, Hg(f) = Hg.) At this point the polariton 
‘A’ is phonon-like and well approximated by the phonon 
mode ‘b’. In the absence of other interactions, the Hamil¬ 
tonian ([3]) accomplishes a perfect quantum state transfer 
between the modes ‘b’ and ‘c’ for an appropriate interac¬ 
tion time. For modes initially in thermal states the result 
can be regarded effectively as a "heat exchange" between 
them. In practice both optomechanical coupling and dis¬ 
sipation are still effective during T 2 , though, so this is 
only approximately correct. However these additional 
effects do not change the dynamics significantly for large 
detunings between the photon and phonon modes and a 
short enough r 2 , as confirmed by numerical simulations 
that account for the full dynamics. 

We now estimate the mean excitation transferred be¬ 
tween the two modes from the state transfer Hamiltonian 
only, noting that In this limit, and for Hg constant during 
T 2 , the dynamics of the two-mode system is governed by 
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the Heisenberg equations of motion 

^ =-i{u}bb + noc); ^ =-iiSc + flob), ( 6 ) 

at at 

which yield readily 

(iV,)(t) = (7V,)(0) + [(iV,)(0)-(iV,)(0)]gsin2(m). (7) 

Here Nj is the number operator for mode ‘j’ and f2 = 
[(wb — 5)^/4 + plays the role of an effective Rabi 

frequency. The maximum exchange of excitation between 
the two modes occurs after the interaction time T 2 = 
7r/(2H). 

The third, compression-like step is an adiabatic change 
that returns the polariton mode ‘A’ to its photon-like 
nature by changing A/ back to A^. As for the previous 
adiabatic step, this must take place in a time T 3 slow 
enough to guarantee adiabaticity, but fast enough that 
thermal relaxation remains negligible. The mechanical 
interaction is off again for the remainder of the cycle. 

The last step is the thermalization of the now photon¬ 
like polariton with its reservoir, effectively at Ri 0, 
over a time T 4 ^ k~^. Repeating the full cycle allows 
one in principle to achieve the ground state cooling of 
mode ‘c’ from room temperature. 

Numerical simulations. This intuitive description of 
the cooling cycle is confirmed by full numerical simula¬ 
tions of the master equation 

+ Hst,p] + nCaip] + 7 {aw + A[p]} , 

( 8 ) 

that includes all three modes. Dissipation is described 
via the super-operators 

Ca[p] = {na + l)[a'^ap - apa^ -I- h.c.) 

+na[aa^ p — a^ pa+ \v.c.), ( 9 ) 

and similarly for and Cc[p\^ where Uj is the average 
thermal occupation of mode j. 

Results of such a simulation are summarized in Fig. [T] 
In this example, the initial phonon number in the mode 
‘c’ to be cooled is (Ac)( 0 ) = Uc = 12, a relatively low 
value chosen for computational convenience. The po¬ 
lariton mode ‘A’ is photon-like, with ria = 0 at optical 
frequencies., but in this example it is given the unrealis¬ 
tically large value {Na){0) ~ Ua = 0.5 to illustrate the 
role of thermal photons. In contrast the polariton mode 
‘B’ is phonon-like, with an initial mean phonon number 
{Nb){^) = 2. It is often the case that the laser-cooled 
mechanical mode is the lowest frequency one, so we have 
taken ujc > oJb- The fact that {Nb){0) < (iVc)(0) can 
then be thought of as resulting from the precooling of 
the phonon mode ‘b’. Note however that except for its 
small non-adiabatic and dissipative coupling to mode ‘A’, 
the polariton mode ‘B’ plays no active role in the oper¬ 
ation of the heat pump, so its initial population is of no 
significant importance. 
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Figure 1. Dynamics of the mean polariton populations (Aa) 
(green dashed line) and (As) (dotted red line) and of the 
population (Ac) of the phonon mode ‘c’ (solid blue line) for 
a few cooling cycles deep in the quantum regime. Here ujb = 
5 = 2 X 10^ K = 40, Ai = -6 X 10^ Af = -6 x 10^ 
g = Do = 2 X 10^, Ti = ra = 4 X 10“^, ra = 8 x 10“®, and 
r 4 = 0 . 1 , with frequencies in units of the phonon decay rate 
7 and times in units of 7 ”^. The initial average populations 
are (Aa)(0) = 0.5, (As)(0) = 2 and (Ac)(0) = 12. The 
four dots on the horizzontal axis mark the end of the first 
’expansion’ step and beginning of ‘heat exchange’, the end of 
heat exchange and beginning of the ‘compression’ step, the 
beginning of thermalization, and the end of the first cooling 
cycle. 

For short times t < ti both (Na) and (Ac) remain es¬ 
sentially constant. During that time A{t) is adiabatically 
changed from A^ to Af, and the polariton and phonon 
mode ‘c’ are uncoupled. The adiabaticity of the ‘expan¬ 
sion’ step is confirmed by the fact that (Na) remains ap¬ 
proximately constant, with modest non-adiabatic transi¬ 
tions between the two polariton modes. 

In the short interval between ti < t < ti + T 2 the in¬ 
teraction between modes ‘b’ and ‘c’ is switched on. Since 
A ft! 5 at A/, we observe an almost perfect "heat trans¬ 
fer" occurring between the the polariton ‘A’ and mode 
‘c’. The mean population of polariton ‘B’ still remains 
essentially unchanged, which confirms that a simple two¬ 
mode model is an excellent approximate description of 
that step. 

The third stroke, of duration T 3 , is a near-adiabatic 
change of A back to A^. Here non-adiabatic effects are 
most apparent in oscillations of {NB){t). This lasts until 
ft ft 0.08 in Fig. [T] Finally, heat dissipation into the 
environment results in the decay of the mean population 
{NA){t) of mode ‘A’, which is now photon-like, to its ther¬ 
mal equilibrium value ria- At the same time, though, the 
phonon mode ‘c’ is also coupled to a thermal reservoir 
and consequently (Ac)(t) slowly increases. Subsequent 
cooling cycles permit to keep it at a value near the quan¬ 
tum ground state, as illustrated by the next two cooling 
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cycles in Fig. [T] (Note that although for the parameters 
of the simulation « 1 optical decay has no detrimen¬ 
tal effect during the ‘expansion’ and ‘compression’ stages 
since it extracts excitations from the system. In addition 
its contribution to non-adiabatic effects is small (^.l 
Cooling limit. To evaluate the cooling limit we first 
consider the "heat exchange" step. For the small nega¬ 
tive detuning A/, the polariton mode ‘A’ consists almost 
entirely of the phonon mode ‘b’. Also, the polariton mode 
‘B’ is far detuned and remains essentially uncoupled, as 
we have seen. Under these conditions we can approxi¬ 
mate the Hamiltonian ([3]) by 

Hst ~ c +c +A), ( 10 ) 

a form that results from carrying out the Bogoliubov 
transformation that diagonalizes iJom, neglecting the 
counter-rotating terms in the resulting interaction be¬ 
tween the polariton ‘A’ and phonon mode ‘c’, and omit¬ 
ting the coupling to polariton ‘B’ altogether. The effec¬ 
tive coupling strength is fig = uflo, with u ^ 1 for small 
negative detunings Ay and g < |A/| 0 is the Bogoli¬ 
ubov transformation coefficient between modes ‘A’ and 
‘b’. We have confirmed numerically that {Nc){t) obtained 
from Eq. cni) and from a full three-mode analysis coin¬ 
cide almost perfectly for our choice of parameters. 

The expression for the mode dynamics obtained from 
the Hamiltonian cni) has the same analytical form as 
Eq. 0, with flo —>■ fig. For the cooling cycle and the 
optimal choice fl't = tt/ 2 , with fl' = [a;yi(Ay) — 5]^/4 -|- 
flg^l^A^ the mean population {Nc,out)j at the end of the 
heat exchange step is related to its value (iVcgn)j; prior 
to that step by 

{Nc,out)j = (1 - + r]{NA,in)j, (H) 

where g = (flg/fl')^ and {NA,in)j is the mean number of 
‘A’ polaritons. If it is properly thermalized while photon¬ 
like at the end of the previous cooling cycle we have 
{NA,in)j = na, independently of j. 

We can determine {Nc,in)j by noting that between heat 
exchange steps the phonon mode is decoupled from the 
polaritons and is only subject to thermalization, so that 

{Nc,in)j =nc + r{{No,out)j-1 - ric), ( 12 ) 

where r = exp[— 7 (ti + T 2 + + T 4 )]. Substituting 

Eq. (I12p into Eq. (Ilip and taking the asymptotic limit 
{Nc,out)j-i = {Nc,out)j gives the cooling limit 

(Aic,out)oo — 7Z T [g‘Oa T (1 ^)(f V)^c\ • 

1 — r(l — g) 

(13) 

The contribution of the phonon thermal noise is fully sup¬ 
pressed for g = 1, which corresponds to the resonance 
condition a;/i(A/) = S. With oJAi^f) ~ Wb, this gives 
LOs = LUc — (jJb, which is precisely the resonance condi¬ 
tion for the substrate modulation frequency to establish 
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Figure 2. Example of two-mode cooling, showing the dynam¬ 
ics of the mean polariton population {Na) (green dotted line) 
and of the populations {Nc) (solid red line) and {Nd) (dashed 
blue line) of the phonon modes ‘c’ and ‘d’ for two cooling cy¬ 
cles deep in the quantum regime. The parametric coupling 
between the modes ‘A’ and ‘c’ occurs at times yt = 0.004 
and 0.0416, and the parametric coupling between ’A’ and ‘d’ 
at times yt = 0.02228 and 0.0604. Same parameters as in 
Fig. 1, except that y = ac/400. The initial average popula¬ 
tions are {Na){0) = 0, {Nc){0) = 10, and {Nd){0) — 7.4. The 
frequencies of modes ‘c’ and ‘d’ are ujc = 1.5iJb and oJd = 2 uJb. 


the state transfer coupling between the phonon modes 
‘b’ and ‘c’. This ideal case yields the fundamental limit 
{Nc^out)oo = na- That is, phonon mode ‘c’ can ideally be 
cooled to the temperature of the electromagnetic field, 
Ta ~ 0 for visible radiation. 

Summary and outlook. In conclusion, we have pro¬ 
posed and analyzed a variation of a heat pump that ex¬ 
ploits the dispersion relation of a polariton fluid to cool 
mechanical modes of arbitrary frequency in an optome¬ 
chanical system. This can also be understood as a form 
of reservoir engineering 0 , [ 2 ^ where changing the na¬ 
ture of a polariton resnlts in its coupling to reservoirs of 
different effective temperatures. Importantly, this cool¬ 
ing scheme can readily be extended to a multiplicity of 
modes by cycling the value of the mechanical coupling 
with the mode ‘b’ from one to the next. Figure 2 shows 
a simulation in which the polariton heat pump is used to 
cool two mechanical modes ‘c’ and ‘d’. 

Although we have considered numerically only the fi¬ 
nal stages of cooling, where the mode starts with only a 
few phonons, the technique can work in principle starting 
from room temperature, although the mechanical mode 
that combines with the optical field to form the polariton 
fluid does need to be colder. A full quantum simulation 
of the process starting from room temperature is beyond 
the capabilities of a numerical simulation, but these first 
stages of cooling could be described classically straight¬ 
forwardly, for instance in the framework of a classical 
Fokker-Planck formalism . 
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